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Abstract. In this paper we examine the rate of convergence of one of the standard algorithms for 
emulating exit probabilities of Brownian motion, the Walk on Spheres (WoS) algorithm. We obtain 
the complete characterization of the rate of convergence of WoS in terms of the local geomnetry of 
a domain. 



1. Introduction 

The harmonic measure on a bounded domain C M*^ at x G can be described as an exit distribu- 
tion of Brownian motion (see [GM04j ) . This measure plays an important role in various problems 
of Probability Theory, Geometric Function Theory, Dynamical Systems, Partial Differential Equa- 
tions, as well as in a vast range of problems of Applied Mathematics. The problem of efficiently 
sampling from harmonic measure is therefore a key problem in Computational Mathematics. 
One of the simplest and most commonly used methods for sampling from harmonic measure is the 
Walk on Spheres (WoS) algorithm. It was first proposed in 1956 by M. Muller in [ Mul56] . Roughly 
speaking, the algorithm consists of replacing the Brownian Motion by a martingale {Xt : t G Z>o}, 
such that Xq = x, and Xt is uniformly distributed on a sphere centered at Xt^i of a radius which is 
a certain proportion of the distance form Xt-i to the boundary dQ (see Section [1.11 for the precise 
definition) . 

It is not hard to see that it takes at most 0(l/e^) steps for the WoS process to reach an e- 
neighborhood of dQ. However, in many situations, this rate of convergence is unsatisfactory. In 
particular, if we wanted to get 2~"-close to the boundary, it would take us a number of steps 
exponential in n. As it turns out, that, depending on the local geometry of the boundary of the 
domain, the rate of convergence is polynomial or even linear in n (i.e. logarithmic in 
Logarithmic rate of convergence of the process Xt to the boundary was established for convex 
domains by M. Motoo in |Mot59] . It was later generalized by GA. Mikhailov in |Mih79] to planar 
domains satisfying any cone condition (i.e. at every point of the boundary there is a cone of certain 
fixed opening in the complement of the domain) , as well for 3-dimensional domains satisfying a cone 
condition with large enough surface angle. See also [EKMSSO] and |Mil95j for additional historical 
background and the use of the algorithm for solving various types of boundary value problems. 
In our earlier work |BB07j . we established polylogarithmic, but not logarithmic, upper bounds on 
the rate of convergence of WoS for planar domains, and for a restricted class of higher-dimensional 
domains. Unfortunately, the techniques of [BB07j do not generalize well to general domains in 
higher dimensions. 

Our present results subsume all prior work on the rate of convergence of the WoS. We introduce 
an easily verified metric condition on the domain which provides tight bounds for the rates of 
convergence. We also show that the condition is tight. 
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1.1. The Walk on Spheres algorithm. Let us now define the WoS. We would like to simulate 
a BM in a given bounded domain Q until it gets e-close to the boundary di}. Of course one could 
simulate it using jumps of size 5 in a random direction on each step, but this would require 0(1/6'^) 
steps. Since we must take 6 = 0(e), this would also mean that the process may take 0(l/e^) steps 
to converge. 

The idea of the WoS algorithm is very simple: we do not care about the path the BM takes, but 
only about the point at which it hits the boundary. Thus if we are currently at a point Xt G 0, and 
we know that 

diXt) :=d{Xt,dn)>r, 

i.e. that Xf is at least r-away from the boundary, then we can just jump r/2 units in a random 
direction from X^ to a point Xf^i. To justify the jump we observe that a BM hitting the boundary 
would have to cross the sphere 

St = {x : \x-Xt\=r/2} 

at some point, and the first crossing location Xt+i is distributed uniformly on the sphere. There is 
nothing special about a jump of d{Xt)/2 and it can be replaced with any (3 d{Xt) where < /3 < 1. 
Let {7t}, t G Z>o be a sequence of i.i.d. random variables each being a vector uniformly distributed 
on the unit sphere in R'^. We could take, for example, 7t = rf/|r^|, where is a normally 
distributed d-dimensional Gaussian variable. Then, schematically, the Walk on Spheres algorithm 
can be presented as follows: 

WalkOnSpheres(Xo, e) 

n := 0; 

while d{Xt) = d{Xt, dn) > e do 

compute rj: a multiplicative estimate on d{Xt) such that /? • d{Xt) < Vf < d{Xt); 

Xt+i := Xt + (n/2) • 7t; 

t:=t + l; 
endwhile 
return Xt 

Thus at each step of the algorithm we jump at least (3/2 and at most 1/2-fraction of the distance to 
the boundary in a random direction. An example of running the WoS algorithm in 2-d is illustrated 
on Figured) 




Figure 1. An illustration of the WoS algorithm for d = 2: one step jump (a), and 
a possible run of the algorithm for several steps (b) 

The proof of the convergence of the algorithm can be found, for example, in |GM04| . Moreover, it 
is not hard to see that the WoS process hits the e-neighborhood of dQ in 0(l/e^) steps. However, 
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in many situations, this rate of convergence is unsatisfactory. In particular, if we wanted to get 
2~"-close to the boundary, it would take us a number of steps exponential in n. As it turns out, in 
many natural situations, the rate of convergence is polynomial or even linear in n (i.e. logarithmic 
in 1/e). The object of the paper is to prove that this is the case, and give precise condition on 
when the faster convergence occurs. 

While an actual implementation of the WoS would involve round-off errors introduced through an 
imperfect simulation, we will ignore those to simplify the presentation as they do not affect any 
of the main results. Thus the problem becomes purely that of analyzing the family of stochastic 
processes {X^} and their convergence speed to d^l. 

1.2. Results. Let Hj^^K) denote the f3- dimensional Hausdorff content of K. 

H0(K)= inf 



A- 



Definition 1. A domain Q (ZW^ is said to be a-thick < a < d if there exists a constant C > 
such that for every x G dO. 

H'^-''{B{x, r)\n)> Cr'^-'', r < 1 

Roughly speaking, a-thick domains have complements of codimension a, which are uniformly large 
at every scale at every boundary point. 

We call the constant c the thickness of the domain $7. It is not hard to see that the property of 
a-thickness is monotone: an a-thick domain is a'-thick for a < a' < d. 
Let us list some examples of a-thick domains. 

(1) All d-dimensional domains are d-thick; 

(2) all bounded d-dimensional domains 0, such that the complement ^1'^ is connected are d — 1- 
thick. 

(3) all convex domains and all domains satisfying cone condition are 0-thick; 

(4) all domains 0, that are bounded by a smooth hypersurface dfl are 0-thick. 

It turns out that the a-thickness of the domain is responsible for the rate of convergence of the 
WoS algorithm. This idea is formulated precisely in our Main Theorem. 

Theorem 2. Let he a hounded a-thick domain in W^. Then the expected rate of convergence of 
the WoS from any x ^VL until termination at distance < e to the houndary is given by the following 
table: 



(1) 





Rate of convergence 


a < 2 


O (log 1/e) 


a = 2 


0(log^ 1/e) 


a > 2 


0((l/e)^-^/") 



The O(-) in the expressions above depends on the dimension d, on a, on the thickness constant C 
from Definition^ and on /3 > from the definition of the WoS. It does not depend directly on il.. 
Moreover, the rates of convergence above are tight. That is, for each a there is a family of a- 
thick domains 17" with some thickness C, such that the rate of convergence with e = 1/n on fi" is 
asymptotically given by the formulas in ([T|) . 

The rate of convergence cannot be better than 0(log 1/e) since at each step of the WoS, the distance 
of Xt to the boundary dO, decreases by at most a multiplicative constant. An intuitive explanation 
to the phase transition phenomenon occurring at a = 2, is that a BM in almost surely "misses" 
sets of co-dimension > 2, while hitting sets of co-dimension < 2 with positive probability. 
It is worth noting that the main result in |BB07j is the special case a = d = 2 of the theorem. 
The following corollaries are implied directly by the Theorem [2j 
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Corollary 3. (1) Since any planar domain is 2-thick, the WoS converges in 0(log^ steps; 

(2) since any planar domain with connected exterior is 1-thick, the WoS converges in 0(log 
steps; 

(3) since any domain in is d-thick, for d > 3 the WoS converges in 0{{\/ e)'^~^/'^) steps; 

(4) since any 2>- dimensional domain with connected exterior is 2-thick, the WoS converges in 
0(log^ 1/e) steps; 

(5) since for any d > 4, any d-dimensional domain with connected exterior is d — 1-thick, the 
WoS converges in 0{{\/ ef'~^/^'^~'^'^) steps; 

(6) since any domain hounded by a smooth hypersurface is 0-thick, the WoS converges in 
0(log 1/e) steps. 

The rest of the paper is organized as fohows. In Section [2] we construct the auxihary boundary 
barrier measures and the energy functions. Using these functions, we prove the upper estimates 
of Theorem [2j More technical estimates on the energy function are done in Section [3l Finally, in 
Section HI we present examples of a-thick domains with the slow rate of convergence of the WoS 
process. 

2. Upper bounds: energy functions 

2.1. Construction of an auxiliary measure. In this section we will construct a family of mea- 
sures near boundary points of an a-thick domain. These measures will be used to construct energy 
functions, which, in turn, play crucial role in the proof of Theorem [2j 

Lemma 4. There exists a constant c = c{a, d, C) such that for any a-thick domain Q. with thickness 
C in and for any x G d^, one can find a Borel measure /x^. which satisfies the following 
conditions: 

(1) sn]iY>{iix) n O = 0, or, equivalently, Hxi^) = 0; 

(2) for any y eM.'^ and r > 0, ^x{B{y,r)) < r^-°; 

(3) for any r <1, /j,x{B{x, r)) > c ■ r'^^"'. 

With a slight abuse of notation, we will also refer to the constant c from this lemma as the thickness 
of the domain. 

Proof. The proof of the Lemma follows the standard reasoning that can be found in, say. Chapter 
II of |Car67j . 

Let us consider the dyadic grid selected so that the point x has coordinates (1/3, 1/3, • • • ,1/3). For 
an integer d-multi-index 7 = (71, ... , 7^), let D^^-y be the cube 

{(xi, ...,Xd : 7n2~'' < Xn < (7n + 1)2"^, n = l,...,d}. 

Let Dk{x) be the unique dyadic cube of the size 2~^ which contains x. Note that x is always at 
distance 2~^/3 from the boundary of Dk{x). 

We will construct inductively the sequence of measures Vn- They will satisfy the following properties: 

(a) suppf„ n $7 = 0. 

(b) VniDka) < H'^''^{Dk,y \n) (or 1 < k <n 

(c) Un{Dk{x)) = H'^-°'{Dk{x) \ n) for 1 < k <n 

Let vi be a delta measure in a point of Di{x) \ Q with the total mass H'^~"{Di{x)). It clearly 
satisfies all of our assumptions. 

Assume now that the measure Vn has already been constructed. The measure I'n+i will be a sum 
of delta-measures on the points outside of fl lying in the cubes from the n -\- 1-st dyadic generation, 
such that i'n+iiDn,'y) = Vn{Dn,'y) for all 7 (so Vn+i will be obtained from Un by re-distributing the 
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latter over the cubes of the (n + l)-st generation). Thus the measure fn+i would automatically 
satisfy the second and the third condition for k <n. 
To construct z/^+i, we use the following rule. 

First, we set Vn+i{Dn+i{x)) = H'^^'^{Dn+i{x)\^). The measure Vn+i clearly satisfies our condition 
(jcj) on Dn+i{x). 

Second, for any other dyadic cubes Dn+i^-y C Dn{x), we assign the mass 

„ /r. \ ffd-afn \ n\ (UniDnix)) - l^n+liDn+lix))) 

[^) l^n+l[-L'n+l,'y) = n \ 



SO that Un{Dn{x)) = Vn+iiDnix)). By sub-additivity of the Hausdorff content, 

D„+i,sCDn{x) 

and hence 

(3) Un+l{Dn+l,^) < H^-^{Dn+l,^ \ Q) 

for D„+i,^ C Dn{x). 

Finally, for any other dyadic cubes from (n + l)-st generation, we set 



(4) l^n+l{Dn+l,y) = i?'*"" (Z^n+l,^ \ 



where Dn is the unique cube from the n-th dyadic generation containing Dn+i^-y- Using the sub- 
additivity of the Hausdorff content, as above, we get the estimate ^ for all cubes of the (n + l)-st 
generation. The construction again satisfies Un+iiDn) = UniDn)- 

Let now u be any weak* limit point of the sequence v is still supported outside of Vi. By the 
second property of the measures Vn, 

(5) v{Du,^) < \ f^) < S'^ H''-"{Dk,^) < 3' (Vd)^ ^-^^"'''^ 

for all k. Using the third property of the measures Vn, the a-thickness of and the fact that 
Dj^{x) contains the ball of the radius /?>, we get that 

(6) y{Dk{x)) > H^-^'iDkix) \ n) > cS-"^ 2-^('^-") 
for any k. 

Every ball can be covered by certain (d-dependent) number of dyadic cubes of comparable size, 
so ([5]) implies that v{B{y,r)) < r'^~°^. Every ball centered at x also contains a dyadic cube of 
comparable (again, d-dependent) size, hence by ([6]), v{B{x,r)) > r'^~^. Now we can set to be 
an appropriately normalized measure v. 

□ 



2.2. Energy Function of optimal growth. The heart of the proof of the upper bounds in 
Theorem [2] is the construction of a subharmonic function with optimal growth at the boundary, 
the Energy Function U on Jl. We will construct U{x) so that it is "small" in the interior of 
and grows to oo as x approaches the boundary d^. The a-thickness of the domain allows us to 
establish that the value of U{Xt) grows in expectation as the WoS progresses. Thus after a certain 
number of steps U{Xt) will be large in expectation which would imply that Xt is close to dfl with 
high probability. 
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The construction of the function is based on the notion of a Riesz potential. For a finite Borel 
measure /x on M"^, and a < d, the a- Riesz potential of the measure /i is defined by 

a — a J \z — x\"- °^ 
For a = the d-Riesz potential is defined by 

[/^(x) = / \og-^dii{z). 
J \z-x\ 

The value Ua{x) = c« is allowed when the integral diverges. 

An important special case is the case of a = 2, the so-called iVewton potential. We will denote U2 
simply by U^. In this case the expression under the integral is harmonic in W^. It is well known 
(e.g. see |Lan72j ) that the function is superharmonic on M*^, and harmonic outside of supp/u. 
More generally, outside of the supp /i, we have the identity 

(7) AU^iy) = - a + 2)(2 - a)U^_,{y). 

It shows that for < a < 2, the function Ua is subharmonic outside of supp fi. 
The following important technical identity, which easily follows from Fubini's Theorem and substi- 
tution, relates the local behavior of the measure /i and the growth of its potential Ua- For a < d, 
we have 

1 r ...o... rKB{y,r)) 



(8) ^"^y^ = d^i ,.{B{y,t'^l^'^--^))dt = ^''^^dr, 
and for a = d, 

(9) C/^(y) = r e-*)) dt = T l^^^^^^lA dr 



Let us now fix an a-thick domain Q. (Z B{Q,1) (ZW^ . Let us consider the set M. of all Borel measures 
II supported inside i?(0, 2) and outside of $7 (i.e. /i(J7) = 0), satisfying the following condition: 

(10) for any y G M*^ and r > 0, /i(S(y, r)) < r'^"" 

Let us now introduce the Energy Function U{y). Recall that U^{y) := U2{y). 

Cin 7-7/„,^ / sup^e^ [/^(y), when q < 2 

^ ' I sup^eA4 U^y), when a > 2. 

Since the set M is weakly* -compact, for every y S there exists a measure maximizing the 
potential in (fTTI) at the point y. 

Let us summarize the properties of U (y) in the following claim. The proof uses the identities ([8]) 
and dH). Recall that d{y) = dist(y,5il). 

Claim 5. Let Q be an a-thick domain. Then 

(1) U{y) is subharmonic in Q. 

2 

(2) For a<2, U{y) < log -— for allyeQ. 

(3) For a>2, Uiy) < -^—diyf''' for all yen. 

a — 2 

Proof. Let a < 2, y £ Q and fi £ Ai. Equations ([8]) and Q imply that 
(12) Uli{y)< C -dt = \og ^ 



d{y) t d{y) 
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Similarly, for a > 2, we will use the harmonic potential U^. Let a > 2, y £ Q and /i E A^. Once 
again, p^ . supp n 17 = 0, and the equations dS]) imply that 



By equations (jl2p and (jl3p . ^7(y) is a supremum of a locally bounded family of subharmonic 
functions. Thus U (y) is subharmonic. 

The second and third statements of the claim follow directly from (jl2p and (jlSp respectively. □ 

Let Xf be the WoS process initiated at some point Xq = y £ Q. Let us define a new process 
Ut = U{Xt), the value of the energy function at the t-th step of the process. Note that because U 
is subharmonic , C/t is a submartingale, that is E[C/j+i|C/j] > Ut- 

For the rest of the section let n = 1/e. Claim [5] immediately implies that a large value of Ut will 
guarantee the closeness to the boundary. More specifically. 

Claim 6. For a<2,ifUt> log2n then d{Xt) < 1/n. 
For a > 2, Ut > {a - 2)n^~'^ implies d{Xt) < 1/n. 

The proof of Theorem [2] relies on finer lower bounds on the function U, which would guarantee the 
optimal rate of boundary convergence. We prove the bounds in the next section. These bounds 
depend heavily on the value of a. We first give a probabilistic proof of the upper bounds in Theorem 
[21 and then prove the finer estimates on U in Section [3l 

2.3. Logarithmic convergence: the case a < 2. In the heart of the proof for this case lies the 
following strong estimate on the behavior of the Riesz potentials near the boundary. 

Lemma 7. For any a < 2 and c > 0, there exist two constants 5 and ij, such that the following 
holds. 

Let $7 be an a-thick domain in M'^ with thickness a. Let y G and x £ be the closest point to 
y. Let fi £ M. 
Then either 

(14) U{z) > C/^(z) + 1 whenever 5/A ■ d{y) <\z- x\ < 5 ■ d{y). 



where x is the point of 50 that is closest to Xt, and p > depends only on (3 and the dimension d. 
Let us fix Xt and take the measure ^ £ M maximizing the value of Ua{Xt). By the preceding 
observation, in the first case in Lemma [71 the subharmonicity of U implies that the expectation of 
Ut-\-k, conditioned on Ut, will increase by some definite constant. 

On the other hand, using the identity ([7]) and the a-thickness of il, one can see that the Laplacian 
of Ua is large near the point Xt in the second case of Lemma [71 Thus, since large Laplacian 
leads to a fast build-up of mean values, we have the above-mentioned increase by a constant after 
the first step. We arrive at the following estimate, which shows that Ut grows at least linearly in 
expectation. 

Lemma 8. There are constants L and k, depending only on c, (3, and a, such that 



(13) 




or 



(15) ii{B{y,2d{y)))>r^d{yY-^ 

The lemma is established in Section [3.11 

Note that after k = 0{\ log5|) steps of the WoS process, 

(16) 5/4 • d{Xt) < \Xt+k — x\ < 6 • d{Xt) with a certain probability p, 



B[{Ut+k-Ut)\Ut]> L. 
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A detailed proof of the lemma can be found in Section 13. 2i 

Lemma [8] implies that E[?7t] > tL/k + Uq. Since d{Xt) > (1 — PYd{Xo), Claim [6] implies that 
Ut < Uo + t\ log(l - + log 2. This implies that Ut > Uq + tL/2k with probabihty at least P, 
where P depends only on /?. This, together with Claim [6] implies the necessary upper bound in the 
case a < 2. 



2.4. Polylogarithmic convergence: the case a = 2. In the case a = 2 the steady linear growth 
of Ut given by Lemma [8] no longer holds. In fact, the only thing that generally holds in this case 
is the submartingale property E[L't+i|C/t] > Ut- We are able to overcome this difficulty by showing 
that the submartingale process {Ut} has a deviation bounded from below by a constant at every 
step. To this end it suffices to show that Ut can grow by some rj with a non-negligible probability. 
We use the following estimate on the energy function (established in Section [3. ip . 

Lemma 9. There exists a constant 6, dependent only on the thickness c, such that the following 
holds. Let Q be a 2-thick domain. Let y € O and x G d^l be the closest point to y. Then 

(17) U{z) > U (y) + 1 whenever \z — x\ < S ■ d{y). 

Since the function U is subharmonic, observation (jl6p implies the following estimate (see Section 
for a proof). 



Lemma 10. Let i7 6e a 2-thick domain in W^. There are constants k and L, depending only on 
the thickness c, the jump ratio (3, and the dimension d, such that 

n{Ut+k-Utf\Ut]> L. 

We can now use Lemma [101 to prove the upper bounds on the rate of convergence for the case 
a = 2. Let us replace the submartingale Ut by a stopped submartingale 



Vt 



t<Tn 

t > T„. 



By the optional stopping time theorem (see |KS9l] ). Vt is also a positive submartingale; Vt < log -. 
This implies, in particular, that 

(18) nVt{Vt+k - Vt)] = B[B[Vt{Vt+k - Vt)\Vt]] > B[E[Vt{Vt - Vt)\Vt]] = 
Lemma [To] implies that 

(19) B[{Vt+k-Vtf]> L-P[Tn>t + k]. 
We are now in a position to establish the upper bounds for a = 2. 

Proof of the upper bound from Theorem [2] for a = 2. Assume first that for some M, 

F[Tn > Mlog^n] > 1/2. 
It means that for all t < Mlog^ n — k, P[T„, > t + k] > 1/2. This implies 

B[V,l,] = B[{iVt+k - Vt) + Vt)^] = E[y/] +E[iVt+k - Vt)^] + 2B[Vt{Vt+k - V)] > B[V,^] + L/2. 

r 1 LAI los^ H 

The last inequality follows from ([HD and ([l9]). Hence E V^^i^ 2„ > rr;^ — . Since V < log ^, 

this leads to a contradiction for large enough M. □ 



2.5. Polynomial convergence: the case a > 2. For the case a > 2, the required converse to 
Claim [6] is relatively simple. 

Lemma 11. For a > 2, and an a-thick domain Q, in M*^ with the thickness c, 

U{y) > K ■ d{yf-^ 
for all y £ 0,. Here the constant K = K{c,a) depends only on c and a. 
The lemma is established in Section [3.11 

The idea of the proof of Theorem [2] in this case is now as follows. When the WoS is far from the 
boundary dQ it makes fairly big steps and when it is close it makes small steps. There are not too 
many big steps because the number of big steps of length > e confined to -6(0, 1) is bounded by 
0(l/e^). On the other hand, there are not too many small steps, because a small step means that 
the WoS is very close to dfl, and should converge before an opportunity to make many more steps. 
More precisely, the number of "big" jumps is bounded by the following Claim. 

Claim 12. Let N{e,T) be the number of the jumps in the WoS process before the time t which are 
bigger then e, i.e. 

N{e,T) = #{t<T\ \Xt - Xt^il > e}. 

Then 

P N{e,T) > 4 < 1/4. 



Proof. Note now that because Xt is a martingale, we have 

T 

(20) 1 > E[X|] - Xl 



T 

■1 = E 

k=l 



E[X?] - E[Xl,]) =Y^E[{Xt - 



i=l 



iXt-Xt^,f]>e'B[N{e,T)] 



t=i 



The last equation implies the statement of the claim, by Tschebyshev inequality. 



□ 




Figure 2. The regions from the proof of the upper bounds for a > 2 



To bound the number of small jumps, we denote hj Rq C 0, the 1/n-neighborhood of 90, and more 
generally, by 

Rk:={x£n : 2^"V™ < di^^ dQ) < 2^ /n} 
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(see Figure [2]). Note that by Lemma [TTl we have 

(21) 2('=-^)(2-a)^a-2 > jj^y^ > ^^k{2-a)^a-2 

for y £ Rk- Using this fact we prove the following. 

Claim 13. Denote by Vk the number of visits of Xt to Rt before the time T when Xt first hits the 
\ / n-neighborhood of the boundary 90, 

Vk = #{t<T : Xt^Rk}. 

Then 

for some constant C2 = C2(c, d, a, (3) and for any M > 1. 

Proof. Suppose that at some point t, Xt G Rk- We estimate from below the probability that this 
is the last time the WoS visits Rk- 

First of all, with some probability p = p{a,f3,c,) > and for some constant r/, Xt+ri G -Rfc+2iogc! 
i.e. the series of first jump brings us much closer to 90. Consider the subharmonic function 

$(y) = (2n)2-"C/(y)-2'=(2-"). 

in O. Then the process ^{Xt+j) is a submartingale. We stop it at time t + t, t > rj, when either 
the WoS terminates or when Xt+r G Rk (i-e. the process gets back to Rk), whichever comes first. 
If Xt+r is 1/n-close to 90 (but not closer than l/2n), then ^{Xt+r) < 1, by (I2T|). If Xt+r G Rk, 
then, again by (I2ip . ^(Xt+r) ^ 0. Another application of (I2ip implies that if y G ^fc+2iogc) then 
$(?/) > 7 for some constant 7. Thus the probability that the WoS terminates at Xt+r (i-e. we 
never visit Rk again) is at least 

P[Xt+r i Rk] > EMXt+r)] > HXt+^) > p-f. 

Thus the probability that the visit Xt to Rk is the last one is at least p ■ 7. The claim now follows 
from an estimate of the probability of having at least Vk returns to Rk , each of them not being the 
last one. □ 

Claims [12] and [13] together imply the upper bounds on the rate of convergence for q > 2. 

Proof of the upper bounds from Theorem [2] for a > 2. By Claim [TS] for any k, we have 
that 

F[vk-s > C2 ■ 2('=-^)("-2) . (3/2 + s/2)] < 1/43/2+^/2 ^ (^/g) . 2-s_ 
Hence, by union bound Vk-s < C2 ■ 2(^-'*)(°-2) . (3/2 + s/2) for ah s > with probability at least 
3/4. Let k be such that 2'^ ~ n^/". Then, with the probability at least 3/4, we have the total 
number of jumps smaller than 2'^/n bounded by 

k k 

(22) '"k-s < X] '^2 • 2('=-^)("^2) . (3/2 + s/2)< 4C2 • 2'=("-2) ^ AC2 ■ n^'^l''. 

s=0 s=0 

If we take = (Ci +8C2)n2~^/'^ steps of the WoS, (p2]) imphes that at least half the steps would be 
of magnitude at least 2^/n ~ n2/'^~^, except with probability < 1/4. Applying the estimate from 
Claim [12] with e = 2^ jn, we see that with the probability at least 3/4, N{2^/n,t) < ^r? ji}^ k, 
4n2~^/". Hence with probability > 1/2 the WoS terminates after 0(?i2~^/") steps. □ 

3. Boundary behavior of the energy function 

In this section we prove the analytical estimates on the behavior of the energy function that have 
been used in Section [2l 
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3.1. Estimating boundary growth. We start with the easiest case a > 2. 

Lemma 1111 (Section 12. 5^ : For a > 2, and an a-thick domain in M*^ with the thickness c, 

\2-a 



for all y G rj. Here the constant K 



U{y) > K ■ d{yy 
K{c, a) depends only on c and a. 



Proof. Let x be the closest to y point at dQ, and let ^ = ;Ua; be the corresponding measure from 
the definition of the a-thick domains. Then, by the identity ^ and since B(x, r) C B{y, r + d{y)), 



(23) Uiy)>U^iy) 



KBiy,r)) 



d(y) 



rd~l 



dr > 



fi{B{x, r — d{x))) 



rd-l 



dr > 



2d{y) 

fi{B{x,r/2)) 

2d(y) 



rd-l 



> c2 



a—d 



2d(y) 



t'-" > K ■ d{y) 



2-a 



□ 



Unfortunately, in the case a < 2 lower bounds of the type established in the proof of Lemma [TT] 
are insufficient, and we will use finer estimates provided by the following construction. 
Let y be a point in $7, x be the point of that is the closest to y, and be a measure in the class 
M. We construct a new measure S A^, which we call the amalgamation of /x at the point y in 
the following way. 




Figure 3. Construction of the amalgamation v = X^fc/^fc 

Let measure be the measure from Lemma H] restricted to B{y,2d{y)), /i2 = = 0, and for 
k > 4, let fik be the measure fj, restricted to the d-dimensional annulus 

Ak = {w : 2'=-id(y) < \w - y\ < 2" diy)} 

scaled by the factor 1 — 7^ := 1 — 2(*^~*'')('^~"). Let us also put 71 = 72 = 73 = 1. We define 

k 

The ingredients of the construction are illustrated on Figure [3l 
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Let us now prove that v G Jv[. Consider any disk B(w,r). Let K be the largest number such that 
B{w,r) intersects Ak- If B[w,r) does not intersect B{y,2d{y)), the measure z/ is no greater than 
fi on B{w,r), and thus v{B{w,r)) < r'^~°'. H K < 3, v{B{w,r)) < iXx{B{w,r)) < r'^"". For ah 
other cases, r > 2^~^d{y), which, by the choice of 77^-, imphes that ^Kf'^~'^ ^ {'^d{y))'^~°^ . Thus 

K 

v{B{w, r)) < fi.iBiy, 2d{y))) + ^,{B{w, r)) - ^ 7kf^{B{w, r) D A^) < 

k=l 

(2d(y))'^-" + (1 - jk)KB{w, r)) < jkt"-'' + (1 - )^''-" = r"-''- 

The second inequahty follows from the fact that the sequence {7^} is non-increasing. We first apply 
the amalgamation construction to the case a = 2. 

Lemma [9] (Section 12. 4^ : There exists a constant 6, dependent only on the thickness c, such that 
the following holds. Let be a 2-thick domain. Let y £ Q and x S dQ be the closest point to y. 
Then 

(24) U{z) > U{y) + 1 whenever \z - x\ < 5 ■ d{y). 



Proof. Since is a compact set, U{y) = U^{y) for some fi £ M.. 

Let fiQ be the restriction of the measure fi to -6(0, 2) \ B{y, 2d{y)). By ([8]) and (fTOl) 

(25) C/^°(y) > C/'^(y)-log2. 
Let V be the amalgamation of /x at y. Next, we will show that 

(26) U''{z)>U^'''{z) -Ci + c- 2^-'^ log \ 



and 

(27) U^'''>{z)>U''°{y) -C2 

whenever Iz — 3;| < Sd{y) for some constants Ci and C2 depending only on d and c. These 
inequalities, together with (j25p . imply the statement of the lemma whenever 6 is sufficiently small 
(namely, when log 1/6 > 2'^-'^{Ci + C2 + 1 + log2)/c). 
To establish ()26p . let us note that for any k we have 

^(^1) + fi{A2) + ■■■ + KAk) = KB{y, 2^d{y))) < {2^d{y)f-\ 
By the Abel summation formula, 

k k 

This implies 

(28) T^E^'^ / I ^-ll^dMw) < Y.lkl^{A,){2^-^d{y)f-'' < 5-4' 

fe>2 ' ' fc>2 
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d-2 



Thus we obtain 

J2Sdiy) r fc>2-^0 ^ 

f'^y^ fi,{B{x,r - 6d{y))) fio{B{z,r)) 1 ^ f dfiojw) 

J25d(y) rd-^ io rd-^ d - 2 ^^^"^ J j,^ \w - z\d-^ 



U%z) > 



> 



which imphes (I26p . 

To obtain ([27]) . note that for any point w £ [y, z], for d > 2 we have, by the estimate (fTU|) 
(29) IV^'^OHI < ^ /|v.^3^| ^/.o(0 =/ '^^0(0 = 

for some constant A depending only on d and c. The same inequahty is derived similarly in the 
case d = 2. This implies that 

f A 

U^'°{z) - W^iy) = / VU'^'iw) ■ dw > -\z - y\—- > -A, 
J[y,z] d{y) 



which is exactly the equation (j27|) . □ 

Another application of the amalgamation construction will establish the lower bounds required in 
the case a < 2. 

Lemma [3 (Section 12. 3^ : For any a < 2 and c > 0, there exist two constants 5 and r], such that 
the following holds. 

Let be an a-thick domain in M*^ with thickness c. Let y G and x G be the closest point to 
y. Let iJ. G M. 
Then either 

(30) Uiz) > U^{z) + 1 whenever 5/ A ■ d{y) <\z-x\<6- d{y). 



or 



(31) M^(y,2%))) >r/d(y) 



Proof. Let rj = 5'^ ""'"^ Assume that ^{B{y,2d{y))) < T]d{yY Let /io be the restriction of ji to 
5(0, 2) \ B{y, 2d{y)), as in the previous lemma. We have, by ([8]), 

(32) UH-iz) > UHiz) - / ^[.IjJ' dr > U^{z) - 2d{y)6. 

J d{z) ' 

On the other hand, the same reasoning as the proof of (j26p above, gives 

(33) U{z) > U^iz) > U^'\z) - Ci + c • 2"-'=' log \ 



for some constant Ci depending only on d, a, and c. Here, as in (j26p . is the amalgamation of ^ 
at y. 

Estimates ()32p and (|33p together imply the statement of the lemma. □ 
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3.2. The boundary drift of the WoS process: a < 2. First we establish that the process Ut 
has the drift toward the boundary in the case a < 2. 

Lemma [8] (Section 12. 3^ : There are constants L and k, depending only on c, (3, and a, such that 

B[{Ut+k-Ut)\Ut]> L. 



Proof. Let us fix Xt. By weak-*-compactness of the set Ai, there exists a measure such that 
U^{Xt) = U{Xt) = Ut. By Lemma [71 either 

(34) U{z) > U^{z) + 1 whenever 5/2 • d{Xt) < \z - x\ < 6 ■ d{Xt). 
where x is the closest to Xt point on dVt, or 

(35) ^^{B{Xu 2d{Xt))) > vdiXtf-'' 

Let us start with the first case. 

For some p > dependent only on d and /3, 

(1 - 0/2)'" d{Xt)/2 < P[\Xt+k - x| < (1 - 0/2)-" d{Xt)] > 
Hence, for sufficiently large k, 

(36) P[6/2 • d{Xt) < \Xt+k -x\<S- d{Xt)] > 

Let us now observe that by subharmonicity of the functions U and Ua, the previous estimate, the 
fact that U > Ua and the assumption ([32]) . 

(37) B[{Ut+k - Ut)\Xt] = E[{U^i{Xt+k) - U^i{Xt))\Xt] + B[U{Xt+k) - U^{Xt+k)\Xt] > 

nU{Xt+k) - U^{Xt+k)\Xt and 6/2 • d{Xt) < \Xt+k - x\ < 6 ■ d(Xt)] > 

Since the value of X^ determines the value of Ut, this establishes the statement of Lemma in the 
first case (with L = p^). 

Now let us consider the second case. By the Green formula, for a C^-smooth function u, 

(38) E[u{Xt+i)\Xt]-u{Xt) = u{y)dS{y)-u{Xt)= r^-M Au{y)dV{y)dr 



lpd(Xt)S'i Jo JB(Xt,r) 

where S"^ is the unit sphere in M'^ with the normalized Lebesgue measure S, and dV is the volume 
element in M"^. 

Note that by (l8|) and ([TD , for | y - | < (3d{Xt) we have 
(39) /\Ul^{y) = {d-a + 2){2-a)U^_M = 

{d-a + 2){2-a)fi[B{Xt,2d{Xt))) H dr > 



l{2+l3)d(Xt) f 

\{2+mxt)r 



for some constant Ci depending only on d, a, and (3. 
So, using ([38]) . applied to n = U^ and ([39]), we get 



(40) nUt+k\Xt] -Ut> B[Ut+i\Xt] -Ut> B[U^'iXt+i)\Xt] - U^^iXt) > 

rPd{Xt) 

r^-\C2{d{Xt))-\'') dr = C2 {{d{Xt) f) {d{Xt))-^ /2 = L 
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for some constants C2 and L depending only on d, a, and (3. Since again the value of Xt determines 
the value of C/j, the Lemma follows. □ 



Let us now turn to the case a = 2. 

Lemma 1101 (Section 12. 4p : Let be a 2-thick domain in W^. There are constants k and L, 
depending only on the thickness c, the jump ratio /3, and the dimension d, such that 

n{Ut+k-Utf\Ut]>L. 

Proof. Fix Xt- By Lemma [9l there exists a constant 5, dependent only on d and c, such that 
(41) U{y) > U{Xt) + 1 whenever \y - x\ < 5 ■ d{Xt). 

This implies that ||f/t+fe — UtW^ > 1 whenever \Xt^k — x\ < 6d{Xt). Note that for some p > 
dependent only on d and f3, 

P[\Xt+k -x\<{l- (5/2)-'' d{Xt)] > p\ 

Hence, for sufficiently large k, P[|Xt4_fc — x| < 5d{Xt)\ > p'^, which, in turn, implies the statement 
of the lemma. 

□ 

4. Lower bounds: examples 

In this section we construct examples of a-thick domains for which the bounds in Theorem [2] are 
tight. The main idea of the construction is as follows. We take a domain A in W^, such as the 
unit ball or a cylinder. We remove a "thin" subset of points C from A to obtain Q = A\C. The 
set C can be thought of as the subset of the grid (7^)"^, for some small 7 > 0. The set C will be 
chosen so that it "separates" the origin from the boundary of A. We set n = 1/e. We choose 7 so 
that the probability of the WoS originated at hitting a 1/n- neighborhood of C before hitting the 
boundary of A is < 1/2 (this means that C is "thin"). Hence, with high probability, the WoS will 
reach dA before terminating. However, in this case the WoS will have to "pass through" the set 
C, where its step magnitudes are bounded by 7. This will, in turn, yield an 0(1/7^) bound on the 
convergence time. The analysis is more intricate in the case when a = 2. In the case when a > 2 
is not an integer, a slight modification to this construction is needed, as will be described below. 

4.1. Proof of the lower bound in the case a > 2. In this section we will give an example of 
a "thin" a-thick domain for which the WoS will likely take r2(n^~^/") steps to converge within 
s = 1/n from the boundary dQ,a- The domain Qa will reside in M*^, where d = \a\ > 3. It is 
easy to see that the examples in higher dimensions d' > d can be constructed from by simply 
multiplying fla by[— 1, l]'' 
The set 

:= (i?(0,l)d_i X [-1,1]) \ 5 

is comprised of a d-dimensional cylinder with a set of points S removed. Here -6(0, l)d-i denotes 
the unit ball in M'^^^. We take A to be the "middle 1/3" shell of the (i-dimensional cylinder: 

A = {zG W^-^ : 1/3 < \z\ < 2/3} x {x G [-1, 1] : 1/3 < |x| < 2/3}. 

Let < 7 ^ 1 be the grid size that will be selected later. We consider the set A^ of gridpoins in 
A. 

A^ = (ci-Lf n A. 
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Let < T] := d — a < 1. Denote by the //-dimensional Cantor set in the interval [0, 1]. It is 
obtained by removing the middle A-fraction of the interval, then removing the middle A-fraction of 
each subinterval etc. For the set to be //-dimensional, we choose A so that 

^ log 2 

log2-log(l-A)' 

In the special case when r/ = 0, we set Co = {0}. We now define the set S: 

S:=A^ + {0} X jG^. 

In other words, S is obtained by attaching a 7-scaled copy of to each gridpoint of A^. This 
completes the definition of the set Qq, = (-^(0, x [—1, 1]) \ S. Each point in dQ^ has an 

r/-dimensional Cantor set in Q'^ attached to it is captured by the following claim. Thus there is a 
universal constant C > 1/16 such that for every 7, the set 0^ is a-thick with the thickness C. 
The following two claims assert that for an appropriately chosen 7, the WoS originated at the 
origin G M'^ and terminated at the 1/n neighborhood of 557^ is likely to hit the boundary of the 
external cylinder (as opposed to the neighborhood of S), and is likely to spend 0(n^~^''") steps 
getting there. 

Claim 14. If 'j > 8n^/"~^ then a WoS originated at and terminated at the 1 / n-neighborhood of 
the boundary dVta will hit the boundary of the cylinder -6(0, Vjd-i x [—1) 1] with probability at least 
3/4. 

Proof. It is not hard to see that we can choose a finite subset P of points in S such that \P\ < 
2j~'^-n^, and for every x such that d{x, S) < 1/n there is ap G P such that |x— p| < 2/n. Consider 
the harmonic function 

(42) ■■= E u_L-2 > 0- 

yep 1^ ^1 

Since the function <I> is harmonic, its application to the WoS process Xt gives a martingale. Hence 
if T is the stopping time of the process, 

E[$(Xt)] = $(Xo) = $(0) < 3'^-^ ■ \P\ < 67-" • n''. 

On the other hand, if d{XT, S) < 1/n, then there is a y G P with \Xt — y\ < 2/n, and 

<1>{Xt) >= l/\y - Xt\''~^ > {n/2f-\ 

Hence the probability of Xt being near S is bounded by 

E[^{Xt)] 67-'^ • n*? 2'^+i7-° 8"n2 
(n/2)rf-2 ^ (n/2)rf-2 ^ n°-2 ^ 4(7n)" ^ ^' ' 

The last inequality follows from the condition on 7. □ 

Claim 15. There is a universal constant 5 > such that for 7 as above, with probability at least 
1/2 the WoS takes at least ^(1/7)^ steps to reach the boundary of the cylinder B{0, l)d-i x [—1, !]■ 

Proof. The proof is done analogously to the proof of Claim [TSl below. □ 

Hence the expected number of steps is at least 

6 / 1 



2 V 8n2/°-i 



17(n2-4/°), 



which completes the proof of the lower bound for Theorem [5] in the case when a > 2. 
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4.2. Proof of the lower bound in the case a = 2. We will now give an example of a two dimen- 
sional domain Q such that the expected convergence time of the WoS to a 0(l/n)-neighborhood of 
dQ is il(log^ n). By taking the d-dimensional domain = Q x for d > 2, we obtain a lower 
bound of r2(log^ n) for 2-thick domains in M'^, proving the lower bound for a = 2 in Theorem [21 
The domain 17 will consist of the unit disc in with O(logn) holes "poked" out of it in a grid 
formation. More specifically, let 7 = 4/ log"^''^ n. We consider the grid T = 7Z x 7Z C M^. We take 
0, to be the unit disc with points from F removed from the "middle third" annulus of the disc. 

(43) n = B{0, 1) \ {{B{0, 2/3) \ B{0, 1/3)) n F) . 

The set il. is illustrated on Fig. [D^a). 

We will show that a WoS originated at the origin Xq = would require an expected time of r2(log^ n) 
to converge. It is immediate to see that the same lower bound holds for any point Xq G B(0, 1/3). 
We first observe the following: 

Claim 16. With probability at least 7/8, a WoS originated at Xq = that runs until d{Xt) < 1/n 
terminates near the unit circle (and not near one of the holes). 

Proof. Let {aj}^^^ = -6(0, 1) \ 17 be the set of holes in Q. Define the harmonic function 

k 

^z) = Y,^og{2/\z-a,\). 

1=1 

It is clear the ^{z) > for all z G -6(0, 1). For any point u in the 1/n-neighborhood of any of the 
holes, <I>(u) > logn. On the other hand, 

^>(0) < A; • log6 < 2/7^ = (logn)/8. 

If Xt is the WoS process with Xq = terminated at time T when d{XT,dQ) < 1/n, then $(Xj) is 
a martingale. Hence, 

(logn)/8 > <^{Xo) = B[^{Xt)] > P[Xt near a hole] • logn. 

Hence the probability that the WoS terminates near a hole is less than 1/8. □ 

For simplicity, we will assume that at every step fo the process the WoS jumps exactly half way to 
the boundary dil. 

To facilitate the analysis we replace the WoS process Xt on with the following process It. It 
evolves in exactly the same fashion as Xt, except when Yt is closer than 1/n to one of the holes 
in 17. In this case, instead of terminating, the process makes a jump of 1/n in a direction selected 
uniformly at random. The process Yt is guaranteed to terminate near the unit circle. We denote 
the termination time by T. Further, we set Yt = Yt for t > T. Note that if the process Xt does 
not terminate near one of the holes, then the process Yj coincides with Xt- Claim [161 implies that 
this happens with probability at least 7/8: 

Claim 17. P[Xt does not coinside with Yt] < 1/8. 

We define two regions A and B, B C A C ft. We take A to be the union of discs with radius 
r = 7/4 around the holes in 17. We take B to be the union of discs with radius r/2 around the 
same holes. The sets 17, A and B are illustrated on Fig. [4l[|a). 

Let time to be the first time with |Yj| > 1/2. Let t' be the first time afterward with either \Yt\ > 2/3 
or ll^l < 1/3. Our goal is to show that with probability at least 3/4, |to — t'\ = 17(log^n). We 
define a subprocess Zt of Yj as follows. Let {si}^^Q be a subsequence of times s between to and t' 
such that Yg ^ A. We set Zi = Yg-. We further define Aj = Zj — An instance of the process 
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Figure 4. An illustration the sets fi, A and B (a), and a possible sequence of jumps 
in the processes {Yt} and {Zt} (b) 



Zi is illustrated on Fig. \Mjo)- Since is a martingale, and Zj is defined by a stopping rule on Yt, 
Zi is also a martingale, and 

(44) E[Ai I Ai,A2,...,A,;_i] = 0. 

In addition, it is not hard to see from the definition of Yt that | Aj| < 4/ log^^^ n for all i. Our first 
claim is that the number k of steps Zi is r2(logn). 

Claim 18. P[k < 10"^ log n] < 1/8. 

Proof. Denote i = 10"^ log n. Then, by 



E[(Zo - Zef] = E[(Ai + A2 + . . . + A,)^] = ^ EfA^] + ^ B[AiAj] = 

j=l l<i<j<^ 

^E[A2]+ E[Ai-E[Aj|Ai]] = ^E[A2] <£. 16/ log n< 1/288. 

i=i i<i<i<<? i=i 

On the other hand, by definition, |Zo — > 1/6, and (Zq — Z^)^ > 1/36. Thus, 

V{k <l\= P[Zi = Zk] < (l/288)/(l/36) = 1/8. 

□ 

Thus the number of steps the process Zt takes is at least 10~'^logn w.p. > 7/8. The process Yt 
consists of the steps of the process Zt plus, in addition, steps the process takes within the region 
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A. We claim that once the process Yt enters the region A, it is expected to spend Q{logn) steps 
there. Moreover, the following holds. 

Claim 19. Let ij > 2. Then there is a 9 > such that whenever Yt £ A, if s > t is the first time, 
conditioned on Yt such that Ys ^ A, then 

(45) P[s - t > 6'log^n] > r//logn, 
for sufficiently large n. 

Proof. Denote the hole in that is closest to Yt by x. Given that Yt £ A, there is some fixed 
probability p > that l^+i € B. In other words, \Yt+i — x\ < r/2 = 7/8. Consider the harmonic 
function 

^{z) = log(r/|x - z\). 

Let t' > t + 1 be the first time such that either Yti ^ A (and thus t' = s), or |yj/ — x\ < n^'Pl^^'^\ 
If Yti ^ A, then ^{Yf) < 0. In the other case, ^iYti) < p\ogn/{4rj). Since t' is a stopping time, 
the optional stopping time theorem applied to the martingale ^{Yt-^-r) combined with the estimate 
^{Yt) > 1/2, gives 

P[\Yt> -x\< > (l/2)/(plogn/(4r/)) > 2r?/(plogn). 

To complete the argument, we claim that assuming \Yt' -x\ < n-P/(5"), it win take the process 
another ri(log^ n) steps to escape A with probability at least 1/2. We consider the process (j)r = 
^{Yti^T-) stopped at time tq when either Vf+ro escapes A, or gets closer than distance 1/n from x. 
0T- is a martingale. Moreover, it is not hard to see that \(po — 4>to \ > plogn/(6r/), and |0j — 0j+i| < 1 
for all i. These two facts imply that 

TO 

E[ro] > Y^Bii^i - <Pi_if] = B[i^ro - 00)'] > iplogn/{6rj)f = / log^ n/(36r?2). 

i=l 

Tschebyshev inequality implies that 9 = /{72rf) satisfies the statement of the claim. □ 

By Claim [TSl we know that except with probability < 1/8 the walk will contain at least r2(logn) 
visits to A. It remains to use Claim [T9l to show that at least one of these stays must be Vt{\o^ n) 
long. Recall that T is the stopping time of the process ly, and k is the number of steps Yt takes 
outside of A. 

Claim 20. Let ai = 10~^ from Claim {TR. There is a constant 02 > such that 

(46) P [fc > ai log n and T < 02 log^ n] < 1/8. 

Proof. For every t such that 1/3 < \Yt\ < 2/3 and Yt ^ A, there is a probability pi > such that 
either If+i G A 01 lt+2 £ A. By Claim [19] we can choose 02 > such that whenever Yf E A, the 
process Yf^r does not escape A for at least a2log^n with probability at least p2 = 6/(aipi log n). 
Hence for each 1/3 < < 2/3 with Yt ^ A, the probability that l^+i or 1^+2 enters A, and stays 
there for at least 02 log^ n steps is at least pi ■p2 = 6/ (ai log n). Since there are at least k = ai log n 
Yts satisfying 1/3 < \Yt\ < 2/3 and Yt ^ A, the probability that for neither one of them does l^+i 
or 1^+2 enter A, and stay there for at least a2 log^ n steps is at most 

(1- 6/(02 log n))'=/2 < (1 -6/(ailogn))"^^°S"/2 ^ g-(6/(ai logn)).(ai logn/2) ^ ^-3 ^ ^/g^ 

□ 

Claims [ITJ [18] and [20] imply the following. 
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Claim 21. Let Xt be the WoS process on the set Q with Xq = 0. Let T' be its termination time. 
Then 

P[T' > azlog^n] > 5/8, 
where 02 > is the constant from Claim [2^ In particular, this implies that E[T'] = f2(log^n). 

Proof. We know that T' > a2log^ n if the following three conditions hold: (CI) the process Xt 
coincides with the process It; (C2) the process Yt makes at least k > ai logn steps outside of A in 
the {z : 1/3 < \z\ < 2/3} annulus; and (C3) the stopping time T of Yt satisfies T > 02 log^ n. In 
fact conditions (CI) and (C3) suffice. We have_P[CT] < 1/8 by Claim [HI P[C2] < 1/8 by Claim 
\IE\ and P[C2 n C3] < 1/8 by Claim [20l Here C denotes the complement of an event C. Hence 

P[CTUC2UC3] < P[CT] +P[C2] +P[C72nC3] < 3/8, 

and 

P[r' > aalog^n] > P[ClnC2nC3] = 1 - P[CT U C2 U C3] > 5/8. 

□ 

Claim [21] gives the lower bound for Theorem [2] in the case a = 2. 
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